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Abstract 

In a super symmetry model with an axino as the lightest supersymmetric particle (LSP) and a 
Bino as the next LSP (NLSP), supersymmetric particle production ends up with including two 
Binos, followed by each Bino's decaying into a photon and an axino. Final states are diphoton 
with large missing energy. We have comprehensively studied the implication of 77 + data 
from the ALEPH, CDF II, ATLAS and CMS experiments. No excess over the standard model 
backgrounds can be explained in this model if the Bino NLSP decays outside the detector. Long 
life time of the Bino is possible because of high Peccei-Quinn symmetry breaking scale f a . The 
ALEPH and CDF II data put a very strong bound on f a for light Bino case with m; ^150 GeV: 
the narrow hadronic axion window around f a ~ 10 6 GeV is completely closed. The recent ATLAS 
and CMS data show very interesting exclusion f a > 10 5 GeV for the Bino mass below 700 GeV. 
This is already stronger than the previous laboratory bounds. 
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Impressive performance of the LHC in its early operation escalates our expectation to 
reveal the secrets of the universe. One of the most profound mysteries is the identity of the 
dark matter (DM). In particle physics, the DM is usually explained by a weakly interacting 
stable particle with a well-motivated symmetry for new physics beyond the standard model 
(SM). If this DM particle is produced at a high energy collider, it escapes detection and 
leaves the signal of missing energy. In order to obtain the missing energy information, we 
need to measure the four-momenta of accompanying SM particles. An isolated photon is a 
good candidate for this role. 

Recently the ATLAS and CMS collaborations have reported the search for diphoton 
events with large missing transverse energy ]fir based on the LHC data at \fs = 7 TeV [1- 
4]. No excess above the SM backgrounds has been found. Together with the previous 
results from LEP2 [5-7] and CDF II [8, 9], these events constrain new physics models which 
predict 77 + $t signal. One good example is the minimal supersymmetric standard model 
(MSSM) with the gauge mediated supersymmetry breaking (GMSB). In this model, the 
light gravitino G is the lightest supersymmetric particle (LSP). A supersymmetric (SUSY) 
particle eventually decays into the next-LSP (NLSP), and the NLSP sequentially decays 
into a gravitino and a SM particle. In most parameter space, the NLSP is the Bino, and the 
Bino decays almost exclusively into a photon and a gravitino. With R parity conservation, 
SUSY particles are always produced in pairs. Therefore all the SUSY final states include two 
photons plus missing energy carried by two gravitinos. In this regard, the experiments of 
Ref. [1-9] provide bounds on the lifetime and mass of the Bino in the GMSB models [10-12]. 

Another interesting new physics model for the diphoton events plus large $t is the SUSY 
model with axino LSP [13]. This model is motivated in order to solve the strong CP problem. 
It introduces an axion superfield with the Peccei-Quinn (PQ) symmetry. In addition to the 
MSSM particle contents, we have an axion and its superpartner axino. This axino can be 
very light with mass about 10 MeV [14, 15]. Then as in the GMSB model, the NLSP Bino 
decays into an axino and a photon. The final states include two photons with missing energy. 
The 77+^T data can have significant implications on this axino-LSP and Bino-NLSP model. 
This is our main goal. 

Let us begin with a brief review on the axion field. The strong CP problem arises from 
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the strong CP odd term of 



where G" is the field strength of a gluon. The absence of neutron electric dipole moment 
leads to extremely small value of 8: \8\ < 0.7 x 10~ n [16]. It requires to be explained by 
some symmetry argument. One way is replacing 9 as a dynamical field 9{x) = a(x)/f a , 
where a(x) is an extremely weakly interacting pseudo-scalar field, called an axion, and f a is 
the axion decay constant. The effective Lagrangian for the axion field is 

C a = -(<V) 2 + -^G a » v G a llv + C^-^F^F^A , (2) 
2 V " ; f a \32n 2 ^ u 7 32vr 2 J ' v ; 

where |C a7 | ~ 1 is a model-dependent parameter. 

After the QCD phase transition, the axion acquires the mass of 
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where z = m u jm^ f n is the pion decay constant, and m n is the pion mass. The axion-photon 
interaction can be rewritten as 

La-yy T Q>Ffj,uF^ i (4) 

where g a ~ n = aC a ~ f / '(2ir /„). 

In the literature, there are two popular models for a very light and weakly interacting 
axion, the hadronic axion model [17, 18] and the DFSZ axion model [19, 20]. Current 
laboratory bound on the axion decay constant is f a > 10 GeV. Much stronger bound comes 
from the astrophysical and cosmological searches such that 10 9 GeV < f a < 10 12 GeV [16, 
21]. Below this cosmological bound, still survives a very narrow but interesting range around 
f a ~ 10 6 GeV, called the hadronic axion window [22]. 

We consider the hadronic axion model with supersymmetry [23]. In our model, the axion 
interaction is described by the following superpotential: 

W = y<S>Q 1 Q 2 , (5) 

where $ = (f) + y/2x9 + F$88 is the axion superfield, and Q± : 2 are SU(2)i singlet heavy quark 
superfields. The vacuum expectation value ((ft) = f a j v^2 is the PQ symmetry breaking scale. 
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The complex field consists of the real-scalar saxion field s and the pseudo-scalar axion field 
a. The axino field is 

' x _y (6) 

The mass of axino is model-dependent. An interesting possibility connected with the DM 
is that the axino can be significantly lighter than other SUSY particles, and becomes the 
LSP [24, 25]. In the no-scale supergravity model, for example, the axino mass is generated 
through one-loop diagrams: 

m« ^ V 2 msvsY, (7) 
167r 

where msusY is the induced SUSY breaking soft mass. If 2/ — 0.1 and msusY — 100 GeV, 
we have m s ~ 10 MeV. If the axino is the LSP, then the NLSP will decay into an axino and 
a SM particle. Hereafter we consider the case where the axino is the LSP and the Bino is 
the NLSP. 

The axino field has the interaction vertices with gg and at one-loop level, mediated by 
two SU(2)l singlet heavy quarks U and D in a simple model. The electromagnetic charges 
of U and D are 2/3 and —1/3 respectively. The effective Lagrangian is 

C~ aSB = " ^fh^B [cos 9 W F^ - sin0 w Zn , (8) 

07T COS Uw J a 

where C a y = 5/3 if the number of U quarks is the same as that of D quarks [21, 23], and 
9w is the electroweak mixing angle. 

At the LHC, axinos are produced as decay products of gluinos and Binos. However the 
decay rate T(g — > ga), which is proportional to 777J/ /f, is extremely suppressed by large value 
°f fa <; 10 4 GeV. Other decay channels of the gluino through strong and/or electroweak 
interactions have much larger decay rates, and thus make Br(^ — > ag) negligible. 

The Bino NLSP is different. Only allowed is its decay into an axino: no matter how small 
its decay rate is, the branching ratio of Br(B — > dj) is almost 100%. The decay width of 
B 07 is given by 
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which leads to the lifetime of Bino as 



lOOGeVN 3 / fJC«Y N 



^ osa& \-^r) [w^v) ns - (10) 

In the axino-LSP and Bino-NLSP scenario, therefore, SUSY particle production leads to 
the final states of diphoton and missing Et, accompanied by SM particles. This process is 
phenomenologically identical to that of the GMSB model. 

As comprehensively summarized in Fig. 1 and 2, we study the implications of the diphoton 
events plus missing energy by the ALEPH, CDF, ATLAS and CMS experiments. All of the 
results are consistent with the SM backgrounds. There are two interpretations for this null 
result. First is that the new physics scale is too high to yield an excess over the backgrounds. 
The new physics cross section is too small. An alternative interpretation is possible in our 
scenario. Binos are produced enough but they decay outside the detector so that we do not 
see 77 + $t signal. This interpretation is possible because of very high scale of f a . In this 
paper, we take the second interpretation. 

We specify our model in more detail. Since main production channels of the Bino pair 
at a hadron collider are the cascade decays from gauginos, gluinos, and squarks, the results 
are sensitive to their mass parameters. For simplicity, we take the following SUSY particle 
mass spectra: 

m^o ~ m~± ~ 1.8 mg, m i R — 1-1 m^, m^ L ~ 2.5 mg, (11) 
m- g ~ 800 GeV, rrig ~ 1.5 TeV. 

The second lightest neutralino X2 an d the lightest charging xt are assumed to be Winos. 
The gaugino and slepton mass relations with the Bino mass are motivated from the GMSB 
SPS8 slope. The gluino mass is from the recent CMS searches for the MSSM signal [26]. 
Note that the slepton masses do not affect the bounds from the CDF, ATLAS, and CMS 
experiments. One possible concern is that our benchmark scenario may have too light 
right-handed selectron. This could be excluded by the OPAL data on di-lepton plus missing 
energy, which provided, at 95% C.L., the exclusion region of the right-handed selectron mass 
and the LSP mass up to mg = 80 GeV [27]. Our condition of ~ 1.1 mg is marginally 
allowed. 

First, the ALEPH data put a very strong bound on f a for light Bino case. The ALEPH 
group searched for the GMSB reaction e + e~ ->• BB ->■ GG77 at y/s — 189 - 209 GeV 
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FIG. 1: Excluded region in the parameter space of (m^, f a ) by ALEPH 77 + CDF 77 + Ifix 
and CDF 7 + ^y+jet data. We set C a y = 5/3 and the benchmark scenario in Eq. (11). 

with a total integrated luminosity of 628 pb" 1 [5, 7]. Crucial event selection criteria is to 
demand photons not originating from the primary vertex of the interaction, which is useful 
for longer Bino lifetime. Our benchmark scenario in Eq. (11) guarantees large enough Bino 
direct production, which is sensitive to the slepton mass. More general slepton masses, but 
still not so heavy compared to the Bino mass, result in similar results. No excess of 77 + $t 
signal over the SM backgrounds can be allowed if the lifetime of the Bino is long enough. 
If the Bino decays within the detection reach, our scenario is not consistent with the null 
result, which excludes small and small f a . In Fig. 1, we show the exclusion region 
denoted by "ALEPH" . The limited cm. energy of the LEP2 covers the Bino mass only up 
to about 94 GeV. In this small Bino mass range, however, the bound on f a is very stringent: 
f a < 10 8 GeV is mostly excluded. The narrow hadronic axion window is completely closed 
for the light Bino case. 

The CDF collaboration also presented their analysis of 77+^+V data, based on 2.6 fb" 1 
total luminosity at y/s = 1.96 TeV [9]. Here X denotes other high-p^ SM particles. This 
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search is motivated by the GMSB model. Significant improvement on the search sensitivity 
was made by a new photon timing system in the EM calorimeter, which measures the arrival 
time of photons [28]. With the tracking device measuring the position and time of the pp 
primary collision, this timing system probes the Bino life time rg up to 2 ns. 

Longer Bino lifetime region in < rg < 40 ns is also covered by CDF experiment through 
7 + jet + lf/T [8]. This final state is motivated by their reference model, the GMSB SPS8 
point. Here the main Bino production channels through gaugino decays are associated with 
prompt taus whose decays can be identified as jets. Photon arrival time is measured by a 
timing system in the ECAL. Final states are a high-i^, isolated, and delayed photon with 
large fa and a high- .Bp jet. Null results with the integrated luminosity of 570pb _1 have 
been reported. 

In our benchmark scenario, gluinos and squarks are too heavy to have sizable production 
cross sections at the Tevatron. The produced Binos are mainly from the cascade decays 
of the gaugino pair production of X2X1 anc ^ XiXi- The inclusive Bino pair production is 
large enough for mg < 150 GeV. If the Bino decays outside the reach of the photon timing 
system, our scenario is allowed. Heavier Bino case, corresponding to heavier neutral and 
charged Winos, leads to too small cross section: the CDF data cannot constrain the model. 
Based on two kinds of CDF data, we exclude the parameter space of (rag, f a ) in Fig. 1. The 
CDF data put also strong bound on f a < 10 7 GeV for 100 GeV < mg < 150 GeV. 

The CMS and ATLAS collaborations [1-3] have also reported their search for diphoton 
and missing transverse energy. With the integrated luminosity of 1 fb^ 1 at y/s = 7 TeV, 
the ATLAS group presented their analysis of 77 + fa final states [2]. The ATLAS search 
is limited for erg < 0.1 mm. At the LHC, multiple collisions from high luminosity make 
it very challenging to measure the timing separation between the primary collision and the 
photon arrival. In our model with s/s = 7 TeV, the gluino pair production becomes also 
important especially when the Bino mass is large. Following our basic interpretation, the 
decay of a Bino outside the reach of erg < 0.1 mm is allowed. This excludes the (rag, f a ) 
parameter space as in Fig. 2. 

The CMS collaboration also presented their preliminary data of 7 + fa + jet event at 
integrated luminosity of 2.1 fb _1 [4]. The CMS group applied the photon conversion impact 
parameter method and searched for long-lived (2 cm < erg < 25 cm) Bino decaying into a 
photon and a gravitino. The transverse impact parameter of a photon to the beam line is 
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FIG. 2: Excluded region in the parameter space of (m^, f a ) by the ATLAS 77 + I^t data and the 
CMS preliminary 7 + ^r+jet data. We assume that the inclusive production cross section of Bino 
is large enough for both ATLAS and CMS. 

measured, which is a robust observable when the true primary vertex is not known because 
of multiple collisions. As in the CDFII delayed photon detection, this method also measures 
the delayed photon from long-lived Bino decay. The CMS preliminary results set the upper 
limits on the Bino inclusive production cross section o < (0.12 — 0.24) pb. In our benchmark 
model, gluino pair production is the main channel with the total cross section about ~ 100 fb. 
Based on these upper limits, we exclude the parameter space of and f a , as in Fig. 2. The 
band structure of the exclusion region is from the limitation of the CMS photon conversion 
impact parameter method. The ATLAS and CMS data extend the search for larger Bino 
mass. The bound on f a is quite significant, which closes a large portion of hadronic axion 
window f a ~ 10 6 GeV especially when the Bino is heavy. 

In conclusion, we obtain the improved bound on the axion decay constant and the Bino 
mass in the axino-LSP and Bino-NLSP model from the ALEPH, CDF II and early LHC 
data on 77 + Ifi? signal. We take a benchmark scenario specifying the SUSY particle masses. 
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Light Bino case with mass below about 150 GeV has a strong constraint from the ALEPH 
and CDF data: the region of f a < 10 7 GeV is completely excluded. The hadronic axion 
window is almost closed in this region. The 1 fb _1 ATLAS data exclud a large portion 
of f a < 10 6 GeV region up to ~ 700 GeV. And 2.1 fb -1 CMS data exclude a region 
around f a ~ 10 7 GeV. This bound is much stronger than the previous laboratory bound 
(/« > 10 4 GeV). 
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